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Abstract 

Modeling wind speed is one of the key element when dealing with 
the production of energy through wind turbines. A good model can be 
used for forecasting, site evaluation, turbines design and many other 
purposes. In this work we are interested in the analysis of the future fi¬ 
nancial cash flows generated by selling the electrical energy produced. 
We apply an indexed semi-Markov model of wind speed that has been 
shown, in previous investigation, to reproduce accurately the statisti¬ 
cal behavior of wind speed. The model is applied to the evaluation of 
financial indicators like the Internal Rate of Return, semi-Elasticity 
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and relative Convexity that are widely used for the assessment of the 
profitability of an investment and for the measurement and analysis of 
interest rate risk. We compare the computation of these indicators for 
real and synthetic data. Moreover, we propose a new indicator that 
can be used to compare the degree of utilization of different power 
plants. 


Keyword: Wind speed, indexed semi-Markov chains, Monte Carlo simu¬ 
lation, Financial Indicators. 

1 Introduction 

Renewable energy is becoming a big player in the world energy production. 
Already in 2009, the world generating capacity through renewable energy was 
of 26% against the 66% coming from fossil fuel [1]. In this scenario the wind 
energy has been growing continuously [2] although the intrinsic volatility of 
wind speed requires accurate evaluation of wind speed data before installing 
wind turbines. Due to its random behavior [3] and the intrinsic difficulty 
of storage [4], many researchers are still working on the quantification of 
available energy and performances estimation of wind projects [5, 6, 7, 8]. 

For this reason the wind speed modeling play an important role in the 
electrical power generation from wind speed. 

Accurate models permit the forecast of future wind speed and energy 
production (see for example [9, 10, 11]. Many scholars have proposed new 
models that can allow the prediction of wind speed, minutes, hours or days 
ahead. Many of these models are based on neural networks, autoregressive 
models, Markov chains [12], hybrid models where the previous mentioned 
models are combined, see [13] and the references therein. Often, these models 
are either focused on specific time scale forecasting, or synthetic time series 
generation. In [14, 15, 16], different semi-Markov models were applied to 
wind speed modeling and it was shown that the semi-Markov framework 
over perform the Markov models and therefore they should be preferred in 
the modeling of wind speed. 

The approach we propose here is based on an indexed semi-Markov chain 
(ISMC) model that was advanced in [17] and applied to the generation of 
synthetic wind speed time series. In [17] we showed that the ISMC model is 
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able to reproduce correctly the statistical behavior of wind speed. The ISMC 
model is a non-parametric model because it does not require any assumption 
on the form of the distribution function of wind speed. Furthermore, it is able 
to forecast wind speed at different time scale without loosing the goodness 
of forecasting which is almost independent from the time horizon [18]. 

The main purpose of this paper is to present a prohtability analysis of 
wind energy production through economical indexes. Indeed, we hrst present 
a comparison between the real energy production and that obtained by im¬ 
plementing the ISMC model on real data. Then we advance a new indicator, 
that can be interpret as a degree of utilization of the wind turbine, and we 
compute it for real and synthetic data. Finally we consider hnancial indica¬ 
tors of the cash flow generated by selling the electrical energy produced in 
the wind farm. The cash flow sequences are unknown because they depend 
on the future random sequences of wind speed that are generated according 
to the ISMC model. The cash flows are subject to the interest rate risk 
which is an important factor that affects the risky operation of producing 
wind energy and the hnancial indicators gives the possibility to manage such 
a risk. 

The paper is divided in this way: hrst, the database and a commercial 
wind turbine is illustrated. Second, the ISMC model is shortly described as 
the theoretical support of the empirical application. Third, a real data appli¬ 
cation is executed involving the proposal of a new index and the computation 
of classical hnancial indicators in our stochastic framework. 


2 Database and commercial wind turbine 

The wind speed data that we employ for the analysis are the same used in 
precedent works [14, 15, 16]. They cover almost hve year with time sampling 
of 10 minutes and the measurement instruments are located at 22 m above 
the ground. Since we need to test the model for the energy production of 
a commercial wind turbine, we transpose our database at the turbine rotor 
altitude using a well known relation, see e.g. [19]: 

(^) “ = i 

where is the wind speed at the height h of the wind turbine hub and 
is the value of the wind speed at the height h^if of the instrument. In our 
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Sate 

Wind speed range m/s 

1 

0 to 3 

2 

3 - 4 

3 

4 - 5 

4 

5 - 6 

5 

6 - 7 

6 

7-8 

7 

8 - 9 

8 

>9 


Table 1: Wind speed discretization 


application h = 50m and hrif = 22m. The symbol Zo denotes a parameter 
that takes into account the morphology of the area near the wind turbine. 
For a region without buildings or trees it varies from 0.01 to 0.001, instead 
for the offshore application it is equal to 0.0001. For our analysis we consider 
a mean value for an onshore application, then we choose Zq = 0.005. 

The ISMC model is based on a discrete state space, for this reason we 
select 8 wind speed values, see Table 1, chosen to cover all the wind speed 
distribution. Table 1 shows the wind speed states with their related wind 
speed ranges. The transformation between wind speed and produced energy 
is made through a commercial 10 kW Aircon HAWT wind turbine. The 
power curve of this turbine is plotted in Figure 1 and its discrete numerical 
values are showed in Table 2. The power curve of a specific wind turbine 
represents the energy production as a function of the wind speed. The con¬ 
tinuous values of the power curve are obtained with linear interpolation of 
the discretized values. 


3 The indexed semi-Markov chain model 

The general formulation of the ISMC model has been developed in references 
[17, 15, 16], here we only discuss it informally. 

Semi-Markov processes have similar idea as those that generate Markov 
processes. The processes are both described by a set of hnite states Vn whose 
transitions are ruled by a transition probability matrix. The semi-Markov 
process differs from the Markov process because the transition times T„ are 


4 













wind speed (m/s) 

Figure 1: Power curve of the 10 kW Aircon wind turbine 


Wind speed [m/s] 

Power [kW] 

1 

0 

2 

0 

2.5 

0.4 

4 

0.9 

5 

2.3 

6 

3.3 

7 

4.9 

8 

6.7 

9 

8.1 

10 

9.3 

11 

9.8 

11.5 to 25 

9.8 


Table 2: Discrete values of the power curve of the 10 kW Aircon wind turbine 


5 





Sate 

jjm j.g] 2 ge m/s 

1 

0 to 2.1 

2 

2.1 - 2.6 

3 

2.6 - 3.4 

4 

3.4- 6 

5 

>6 


Table 3: f/™ processes discretization 


generated according to random variables. Indeed, the time between transi¬ 
tions r„+i — Tn is random and may be modeled by means of any type of 
distribntion fnnctions. In stndies concerning wind speed modeling the states 
Vn indicates discretized wind speed at the nth transition and T„ the time in 
which the nth change of wind speed occnrs. 

In order to better represent the statistical characteristics of wind speed, 
the idea of an ISMC model was advanced in the held of wind speed, see 
[17]. The novelty, with respect to the semi-Markov case, consists in the 
introdnction of a third random variable dehned as follow: 


r) 


k=Q 


T 


n—k 


Tn-l-k 


T - T 

-L ri -1-7 


n—l—m 


( 2 ) 


This variable can be interpreted as a moving average of order m -|-1 execnted 
on the series of the past wind speed valnes (vn-i-k) with weights given by the 
fractions of sojonrn times in that wind speed (Tn-k — Tn-i-k) with respect 
to the interval time on which the average is execnted (T„ — Tn-i-m)- 

The ISMC model considers that the probability of changes in wind speed 
do depend also on this new variable. In a simple semi-Markov model it 
depends only on the present wind speed and on the transition time. The 
ISMC model can be considered as an m order semi-Markov chain where the 
dependence is given only by the averaged last m states. Also the process 
f/™ has been discretized. Table 3 shows the states of the process and their 
valnes. 

The parameter m must be optimized as a function of the specihc database. 
The optimization is made by hnding the value of m that realizes the minimum 
of the root mean square error (RMSE) between the autocorrelation functions 
(ACT) of real and simulated data, see [17]. In our analysis m = 7. The 
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ISMC model revealed to be particularly efficient in reproducing together the 
probability density function of wind speed and the autocorrelation function, 
see [17]. 

The one step transition probability matrix can be evaluated by consid¬ 
ering the counting transition between the three random variables considered 
before. Then, the probability Pij{t,u) represents the transition probability 
from the actual wind speed state i, to the wind speed state j, given that the 
sojourn time spent in the state i is equal to t and the value of the process 
77™ is u. These probabilities can be computed as: 




^i,j (7) n) 
j 


(3) 


where riij [t, u) is the total number of transitions observed in the database 
from state i to state j in next period having a sojourn time spent in the wind 
speed i equal to t and the value of the index process equal to u. 


4 Economic application 

Mathematical models of wind speed can be used to assess the economic valid¬ 
ity of an investment in a wind farm. In this section, we show how the ISMC 
model can be used to evaluate useful economic indicators like the Internal 
Rate of Return, the semi-Elasticity and the relative Convexity of the cash 
flow generated by the investment on a wind farm in a specihc site. We con¬ 
duct two separate analysis: in the hrst one we consider an investment for a 
period of 15 years, which corresponds to the period in which the government 
incentives on the energy price in Italy are still present. In the second case we 
consider an investment without the government incentives and in this case, 
being the price of the energy lower than in the previous case, the period of 
the investment is settled to be greater and equal to 30 years. 

Since the database of wind speed is composed of only 5 years of wind 
speed measures, we apply a bootstrap procedure in order to have the suf- 
hcient number of data for the analysis. The database has been adapted to 
the blade height by using equation 1. To compare real data with the ISMC 
model we generate, through Monte Carlo simulation, 1000 possible scenarios 
of wind speed time series of length 15 years for the hrst case and 30 years 
for the second case. The sampling period is 10 minutes. 
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Real 


ISMC 




Energy [kWh] Energy [kWh] 

Figure 2: Histogram of the electrical energy produced in 10 minutes. Com¬ 
parison between real and simulated data. 

4.1 Energy production 

According to the power curve described in Figure 1 and Table 2 wind speed 
has been transformed into power and then in produced energy in a specihc 
time window. In Figure 2 it is shown the comparison between the distribution 
of the produced energy in a 10 minutes interval for real and synthetic data. 
Although the transformation from wind speed into energy is nonlinear, the 
ISMC model is able to reproduce correctly also a real case of energy produc¬ 
tion. In fact the values obtained in the two cases are almost identical. The 
particular form of these distributions is due to the existence of a cut-in wind 
speed, under which no energy is produced, and a part of the power curve 
that remains with constant value, from 10 m/s until 32 m/s, in which the 
wind turbine produce energy at its rated power. This explain the existence 
of the two modes. 

In table 4 we give some statistics of the annual energy distribution for 
real and synthetic data. We included the results for the JarqueBera test 
(JB) which is a goodness-of-£t test of whether sample data have the skewness 
and kurtosis matching a normal distribution. The test rejects the normality 
hypothesis for the annual energy. This is mainly due to the absence of a 
Gaussian tail of the distribution on its right side (see the distribution figure 












Annual Energy 


Mean 

STD 

Skewness 

Kurtosis 

JB 

jB stat 

p-value 

Real 

147.4 

42.90 

-0.331 

2.70 

Rejected 

29 

0.001 

Synthetic 

150.5 

37.82 

-0.283 

2.65 

Rejected 

24 

0.001 


Table 4: Mean value, standard deviation, Skewness, Kurtosis and JB test for 
real and synthetic energy produced 



Figure 3: Synthetic annual energy distribution as a function of real annual 
energy distribution. 
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in [20]). 

In Figure 3, instead, we show the synthetic annual energy distribution 
against the real annual energy distribution. Dots close to the straight line 
demonstrate the appropriateness of the ISMC model. The coefficient of de¬ 
termination for this set of data is = 0.89. 

As completion of this section, we show another comparison between real 
and synthetic energy produced. Following a similar approach of [21], we 
propose a new indicator that quantify the utilization degree of a wind farm. 
We dehne the satisfying power demand percentage (SPDP) as the part of time 
in which the wind turbine is satisfying the electrical power demand related 
to the specihc considered period. Formally it can be dehned as follow: 


SPDP 


TjPg -Pr>0) 
T 


(4) 


where the numerator represents the period in which the electrical power 
available {Pg) satishes the electrical power required (P^) and T is the number 
of unit of time where the inequality Pg — Pr >0 is satished in the considered 
time horizon. To compute this indicator we evaluate the power available 
(practically the power available on the wind turbine) and we compare it 
with a scaled electrical power demand. We use the data of [21] and scaling 
its maximum value. Particularly in [21] the greatest value (during the day) of 
the power demand is 14.8 MW. We scaled all the curve to have the maximum 
value equal to the rated power of the chosen commercial wind turbine, then 
9.8 kW. The result of the scaling procedure is showed in Figure 4. This 
represents the power demand in a day, we extend it for the considered period 
(hve years) without considering any seasonality. 

In Figure 5 we show the comparison of the SPDP evaluated on real and 
synthetic data for a period of hve years. Moreover we compute it by varying 
the over-dimensioning of the plant. In other words, we hx the maximum 
power required during the day at 9.8 kW and we suppose to oversize the 
rated power of the wind turbine from 0 to 100 % of its value. For example if 
we have a wind turbine with a rated power of 10 kW and we want to improve 
the plant, we are interested in knowing how the SPDP will change. If we 
assume a 40% of over-dimension (then we will have a turbine of 10-1-4 kW), 
the SPDP will be 13%, with a percentage variation of about 10%. 

It is interesting to note that the variation of the SPDP is very important, 
it increases almost of 50% with a variation of 100% of the overdimensioning. 
This behavior is well kept by the model. 
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Hours [h] 

Figure 4: Power demand scaled on the commercial wind turbine rated power 



Overdimensioning [%] 


Figure 5: Variation of the SPDP for a 5 years time period as a function of 
the rated power overdimensioning 
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4.2 Financial indicators 

In this subsection we will compare important financial indicators computed 
on real time series, generated through the bootstrapping procedure, and sim¬ 
ulated data, generated through the ISMC model by means of Monte Carlo 
simulations.We consider two different cases. In the first case we compute 
the cash flow obtained by selling the energy produced by the wind turbine 
taking into account the Italian government incentives on the selling price. In 
Italy for a wind farm with a rated power under 1 MW and until 15 years 
from the setup of the plant, the price of the energy is fixed, at the time of 
writing, to 0.291 €/kWh. In this case we consider a cash flow of 15 years. 
Just to give a quantitative example, the second case takes into account the 
selling of energy at the present market price in Italy, without the government 
incentives, which is 0.027 €,/kWh. In this case the considered period of the 
cash flow is extended to 30 years. 

For each case we generate 1000 possible scenarios and we evaluate the 
financial indicators for both real and synthetic data. The trajectories are 
transformed into power through the power curve of the wind turbine chosen, 
then in energy by considering the time between two successive steps and after 
in euro by multiplying for the specific price considered. The annual cash flow 
are built and at time 0 the cost of the initial investment is considered. Com¬ 
monly for wind turbines with a rated power under 100 kW it is considered a 
price of 3000 € for each kW (this value comes from some quotes requested 
to sector companies). This price takes into account the cost of the turbine, 
the setup cost and also the annual maintenance cost. Given that our wind 
turbine a 10 kW of rated power, we choose the initial investment of 30000 
€. 

The considered indicators are the Semi-Elasticity and the Relative Con¬ 
vexity. We show the histograms of the financial indicators evaluated for all 
the 1000 scenarios, for real and simulated data, and for the two cases with 
and without the government incentives. 

4.2.1 Indexes for the measurement of the interest rate risk 

In this subsection we evaluate some measures of the interest rate risk that 
are commonly used in financial analysis and we apply them to the investment 
in a wind farm. As we have seen, the production of energy generates a cash 
flow sequence. The value of this stream may change in time due to variations 
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Semi-Elasticity 



Mean 

STD 

Skewness 

Kurtosis 

JB 

jB stat 

p-value 

Real 15 

-6.79 

0.0024 

0.0758 

2.92 

Accepted 

1.21 

0.50 

Synthetic 15 

-6.77 

0.0246 

-0.0658 

3.01 

Accepted 

0.72 

0.50 

Real 30 

-7.03 

0.0022 

-0.0183 

3.03 

Accepted 

0.11 

0.50 

Synthetic 30 

-7.05 

0.0118 

-0.0110 

3.10 

Accepted 

0.43 

0.50 


Table 5: Statistics of semi-Elasticity distributions evaluated for real and 
synthetic data and for the two cases: with government incentives (15 years) 
and without (30 years). 


of the interest rate. The wind energy producer, as any other hrm, may suffer 
the risk of changes in the interest rate and then it is necessary to control 
such kind of uncertainty. 

One of the most important index is the Duration (D) that was introduced 
by Mcaulay (1938). It is defined as a weighted average of time to payment 
of the cash flow sequence. From the Duration it is possible to recover the 
semi-elasticity (S'e) that is a relative measure of the variation of the value of 
the cash flow stream. In formula we have: 

. ^_ 1 E:=iS-CF,-{l + r)-^ ^ D 

1 + r ■ (1 +l + r' 


In the application the interest rate r was equal to 3% and the results on the 
semi-elasticity are summarized in Figure 6. In Table 5 instead, we give some 
numeric values of the statistics of the distributions. 

Another important index is the Relative Convexity (RC) that expresses 
the convexity of the discounted value of the cash flow stream in unit of value 
variation. In formula we have: 

Se E:.iCF, ■{l+r)-‘ Se 

and C denotes the convexity of the cash flow. The behavior of the duration 
and of the convexity were computed, by the same authors, in [20]. 

In Figure 7 is plotted the Relative Convexity in the same cases of the 
previous indicators and in Table 6 there are its statistic values. 

As it is possible to note, for each pair of indicators and for each specihc 
case of the period of investment, the mean values of the indicators are almost 
the same for real and simulated data. 
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Bootstrap I SMC 




(a) 


Bootstrap ISMC 



Figure 6: Histograms of the Semi-elasticity for real and simulated data, (a) 
Case with government incentives and cash flow of 15 years, (b) Case without 
government incentives and cash flow of 30 years. 
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Bootstrap I SMC 



(a) 


Bootstrap 


ISMC 


0.12 


O 

c 

CD 

D 

O' 

CD 


0.08 


0.04 



- 17.98 - 17.97 - 17.96 

Relative Convexity 



(b) 


Figure 7; Histograms of the relative Convexity for real and simulated data, 
(a) Case with government incentives and cash flow of 15 years, (b) Case 
without government incentives and cash flow of 30 years. 
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Relative Convexity 



Mean 

STD 

Skewness 

Kurtosis 

JB 

jB stat 

p-value 

Real 15 

-10.52 

0.0025 

0.0768 

3.07 

Accepted 

1.2277 

0.50 

Synthetic 15 

-10.40 

0.0677 

0.0226 

3.24 

Accepted 

2.6253 

0.50 

Real 30 

-17.97 

0.0035 

-0.0024 

2.99 

Accepted 

0.0013 

0.50 

Synthetic 30 

-17.98 

0.0176 

0.0050 

3.16 

Accepted 

1.1506 

0.50 


Table 6: Statistics of Relative Convexity distributions evaluated for real and 
synthetic data and for the two cases: with government incentives (15 years) 
and without (30 years). 


5 Conclusion 

In this paper we apply an indexed semi-Markov model (ISMC) to compute 
different financial indicators like the Internal Rate of Return, the Semi- 
Elasticity and the Relative Convexity. To do this we use a real scaled wind 
speed database and a commercial wind turbine in order to compute the real 
energy produced. To conduce an economic analysis on a medium therm in¬ 
vestment on a wind turbine, we apply a bootstrap procedure. Moreover, 
we propose a new indicator that can quantify the degree of utilization of a 
power plant. By means of Monte Carlo simulation we show that the ISMC 
can reproduce almost exactly all the indicators calculated on real data. 
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